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Abstract 

We have developed a sophisticated model of the radio emission from extensive air 
showers in the scheme of coherent geosynchrotron radiation, providing a theoreti- 
cal foundation for the interpretation of experimental data from current and future 
experiments. Having verified the model through comparison of analytic calcula- 
tions, Monte Carlo simulations and historical experimental data, we now present 
the results of extensive simulations performed with our Monte Carlo code. Impor- 
tant results are the absence of significant asymmetries in the total field strength 
emission pattern, the spectral dependence of the radiation, the polarization charac- 
teristics of the emission (allowing an unambiguous test of the geomagnetic emission 
mechanism), and the dependence of the radio emission on important air shower 
and observer parameters such as the shower zenith angle, the primary particle en- 
ergy, the depth of the shower maximum and the observer position. An analytic 
parametrization incorporating the aforementioned dependences summarizes our re- 
sults in a particularly useful way. 
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1 Introduction 



Radio emission from cosmic ray air showers, initially predicted by Askaryan 
1^ and experimentally discovered in 1965 by Jelley et al. j^, offers the 
opportunity to carry out cosmic ray and neutrino research with radio mea- 
surements in the frequency range from a few to a few hundred MHz. The 
radio technique yields data very much complementary to those collected with 
ground-based particle detector arrays and air shower fluorescence detectors 
and incorporates the advantages of a 100 % duty cycle, moderate cost per 
antenna, and the possibility of observation even in populated areas j^. 

After the flurry of activity in the 1960's and 1970's, research in this field 
virtually ceased completely in the late 1970's due to experimental problems, 
difficulties in the interpretation of the experimental data and the success of 
alternative techniques. With the advent of new fully digital radio interferom- 
eters such as LOFAR^ and the ideas put forward by Falcke & Gorham in 
0], the radio technique, however, currently is experiencing its renaissance. Al- 
though the effect has been known for almost 40 years, the knowledge of the 
radio emission properties and their dependence on the underlying air shower 
parameters has been rather poor until recently [see, e.g., the discussion in|^. 
To address this issue and to build the necessary foundation for the applica- 
tion of radio techniques to cosmic ray research, a number of experimental and 
theoretical efforts have started in the last few years. 



On the experimental side, the LOPES project p is developing and testing the 
design of a LOFAR prototype station for the measurement of cosmic ray air 
showers. It will provide experimental data that should finally pin down the 
absolute strength as well as the polarization properties of the radio emission 
very soon. Another ongoing effort is the CODALEMA experiment in France 

0. 



On the theoretical side, there have been calculations presented in p] and 
Furthermore, we have developed our own analytical model of the radio emis- 
sion from vertical cosmic ray air showers in the scheme of coherent geosyn- 
chrotron radiation in j^. This model allowed us to gain a solid understanding 
of the important coherence effects that shape the radio emission. As a second 
step in our modeling efforts, we designed and implemented Monte Carlo simu- 
lations of the radio emission from cosmic ray air showers in the same scenario, 
yet with much higher precision and for arbitrary geometries . The simula- 
tions are based on analytic parametrizations of the air shower characteristics 
and constitute a precursor to our upcoming full-fledged Monte Carlo simula- 
tions of radio emission from extensive air showers based on precise air shower 
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modeling with CORSIKA [nj. 



In this article, we present the results inferred so far from the simulations per- 
formed with our code. After a short description of the underlying simulation 
parameters in section 2, we describe important characteristics of the radio 
emission in general such as the radial dependence, the spectral dependence, 
the curvature of the radio front and the polarization characteristics of the 
radiation (which play an important role in experimentally verifying the dom- 
inant emission mechanism) in section 3. We analyze the dependence of these 
characteristics on the associated air shower parameters such as the shower ge- 
ometry (zenith and azimuth angle), the primary particle energy, the depth of 
the shower maximum, and the magnetic field in a qualitative way in section 4. 
Afterwards, we parametrize the emission's dependence on the various observer 
and shower parameters in a number of individual formulas for our reference 
shower (section 5) before generalizing the parametrizations to arbitrary shower 
geometries and piecing together an overall parametrization incorporating all 
dependences in section 6. We discuss our results in section 7 and conclude the 
article in section 8. 



2 Simulation parameters 

All simulations presented here are performed with the Monte Carlo code de- 
scribed in [10]. In this section, we specify the simulation strategy as well as 
the parameters and configuration options that are used throughout this work. 

Obviously, we cannot perform a true rz-dimensional analysis of the parame- 
ter space in question. We therefore choose a vertical 10^^ eV air shower as a 
reference and change only one of the shower parameters at a time to analyze 
its effect on the radio emission. This postulates that the effects introduced by 
changes of the different parameters are well-separable. We pay special atten- 
tion in cases where this is obviously not true (e.g., primary particle energy 
and depth of the shower maximum). 

Our reference air shower is calculated with a primary particle energy of 10^'' eV, 
developing to its maximum at an atmospheric depth of 631 g cm^^ as orig- 
inally adopted in j^. This corresponds to a distance of ~ 4 km above the 
ground in the case of a vertical air shower. 

The following settings are kept throughout all simulations if not explicitly 
stated otherwise (definitions of the terms in quotation marks are given in jlQ|): 
The particle track lengths are distributed following an exponential probabil- 
ity distribution with mean track length of 36.7 g cm~^ for both electrons and 
positrons. The particle energies are set to follow a broken power-law distri- 
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bution peaking at 7 = 60 as described in j5|. The magnetic field is chosen 
with a strength of 0.5 Gauss and an inchnation of 70°, which approximately 
corresponds to the configuration present in central Europe. Calculations are 
done on a "simple grid" of 1 ns resolution with "smart trajectory sampling" 
enabled. "Automatic ground-bin inactivation" is used with a precision goal of 
0.25 % in four consecutive blocks of 10000 particles each up to a maximum of 
25000000 particles. A total of 800 bins (32 in azimuth; 25 in radius, up to a 
distance of 1000 m) is calculated in each simulation. 



3 General characteristics 

First, we present the general characteristics of the radio emission from a pro- 
totypical 10^^ eV vertical air shower, which we take as the reference shower in 
this work. 

3.1 Spectral dependence 

In Fig. 1 we present the spectral dependence of the emission from a 10^^ eV ver- 
tical air shower at various distances from the shower center. The spectra show 
a steep decline towards higher frequencies due to the coherence diminishing as 
the wavelength becomes shorter and thus comparable to the scales present in 
the shower "pancake" . The field strength reaches a first interference minimum 
at a distance-dependent frequency. Afterwards, we see a quickly alternating 
series of maxima and minima that are insufficiently sampled in this calculation 
and therefore give rise to the unphysical-seeming features at high frequencies. 
A realistic modeling of the emission in this incoherent regime would need a 
more detailed air shower model taking into account the inhomogeneities that 
are known to be present in the shower cascade. This cannot be achieved with 
the currently used analytic parametrizations of air shower properties, but will 
be accomplished once our code is interfaced to the air shower simulation code 
CORSIKA. 

The larger the distance from the shower center, the steeper the spectral de- 
pendence of the emission. In other words, coherence is much better up to 
large distances at lower frequencies compared with higher frequencies: while 
the emission is coherent to large distances of > 500 m for the 10 MHz fre- 
quency component, it already becomes incoherent at ~300 m for the 55 MHz 
frequency component. 

Figure 2 shows the spectra of the same 10^^ eV vertical air shower plotted in 
a. v Sy diagram, illustrating that most of the power is emitted at frequencies 



4 



around 20 MHz to 30 MHz. 



These two effects strongly point to low frequencies as the most promising 
regime for observation of cosmic ray air showers with radio techniques. 

To demonstrate the differences for emission at low and high frequencies, we 
compare some of the results presented in the following sections for the two 
prototypical frequencies of 10 MHz (good coherence up to large distances as 
desirable for experimental measurements) and 55 MHz (frequency band used 
in the historical works and LOPES, cf. [Sj? but coherence only up to medium 
distances) . 
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3.2 Radial dependence and emission pattern 



Fig. 3 shows the 10 MHz component of the electric field strength in the indi- 
vidual linear polarization directions "north-south" , "east- west" and "vertical" . 
The total field strength pattern is remarkably symmetric in spite of the intrin- 
sic asymmetry of the geomagnetic emission mechanism. A more quantitative 
view of the radial dependence of the emission is depicted in Fig. 4. 

Please note that we can equivalently use the east- west polarization component 
or the total field strength in many of the following analyses as there is no fiux 
in the north-south (let alone vertical) polarization component along the north- 
south direction from the shower center for air showers coming from the south. 
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Fig. 3. Contour plots of the 10 MHz field strength for emission from a 10^'' eV 
vertical air shower. From left to right: total field strength, north-south polariza- 
tion component, east-west polarization component. The vertical polarization com- 
ponent (not shown here) does not contain any significant flux. Contour levels are 
0.25 /iV m^^ MHz"^ apart. 
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3.3 Wavefront curvature 



The radio wavefront arriving at the ground is not planar. As demonstrated in 
Fig. 5, the pulses systematically lag behind at larger distances from the shower 
center. The curvature of the wavefront can be approximated by a spherical 
surface with a given radius. At distances beyond a few hundred meters, this 
approximation, however, breaks down. (Additionally, the curvature radius de- 
rived from the timestamps of the maximum filtered pulse amplitudes depends 
on the specific filter used.) The scatter seen in the plot is not of statistical na- 
ture, but rather represents the slight time-shift in the pulses' peak amplitude 
function of azimuth angle. 

The curvature of the radio wavefront plays an important role for the beam- 
forming performed in digital radio interferometers and seems to be confirmed 
by LOPES measurements (Horneffer, private communication). 
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Fig. 5. Radio wavefront curvature given by c times the time-lag of the east-west 
polarization raw pulses emitted by a vertical 10^^ eV air shower. Lines from top to 
bottom: curvatures given by spherical surfaces with 5700 m, 6700 m and 7700 m 
radius. 

3.Ji- Linear polarization 

The radio emission generated by the geosynchrotron mechanism is intrinsically 
linearly polarized to a very high degree. Figure 6 shows the raw (unfiltered) 
pulses arriving at a distance of 200 m to the north-west from the center of a 
10^'^ eV vertical air shower. The north-south and east-west polarization com- 
ponents are of similar strength and arrive almost synchronously, as expected 
for a linearly polarized pulse. (There is, however, also a small degree of cir- 
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cular polarization since the two peaks arrive not completely in coincidence.) 
The vertical polarization-component is negligible. 

Figure 7 shows the same data (neglecting the vertical component) visualized 
as a scatter plot. For each time-step of the simulated pulse, a point specifying 
the north-south versus east-west field strength component is drawn. In other 
words, the series of points directly illustrates the evolution of the (projected) 
electric field vector. The very narrow "loop" performed by the vector in the 
upper-left quadrant of the diagram demonstrates that the emission is indeed 
linearly polarized to a very high degree, even at the already moderate distances 
presented here. (In case of perfect linear polarization, the series of points 
would all lie on a straight line, whereas for perfect circular polarization, the 
polarization vector would follow a full circle around the origin.) 

In the center regions where the emission is strongest, the radiation is almost 
perfectly linearly polarized. In these regions, the polarization vector points 
in the direction perpendicular to the air shower and magnetic field axes, as 
predicted in [3], cf. section 4.2. 
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Fig. 6. Raw (unfiltered) pulses in the individual linear polarization components at 
200 m distance to the north-west from the center of a 10^'' eV vertical air shower. 
Solid: east-west component, dashed: north-south component, dotted: vertical com- 
ponent. 



4 Qualitative dependence on shower parameters 

In the following subsections, we present a number of dependences of the radio 
emission on specific air shower parameters in a qualitative way. 
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Fig. 7. Scatter plot of the north-south and east-west polarization components shown 
in Fig. 6. The emission is linearly polarized to a high degree. 

4-1 Shower zenith angle 

An interesting question is that of the radio emission's dependence on the air 
shower geometry. Fig. 8 shows the radial dependence of the 10 MHz frequency 
component for air showers coming from the south with different zenith angles. 
It is well visible that the radial dependence in the north (i.e., shower axis) 
direction becomes much flatter with increasing zenith angle. 

A broadening of the emission pattern in the shower axis direction could be 
intuitively expected from projection effects occurring when the air shower is 
inclined. One can remove these projection effects by changing the coordinate 
system from the ground-based "distance to the shower center" to the shower- 
based "(perpendicular) distance to the shower axis". (The electric field vector, 
however, is still denoted with the ground-based north-south, east-west and 
vertical components which thus do not change in strength. This method of 
"back-projection" is not the same as an inclination of the full ground-plane.) 
Fig. 9 shows the back-projected radial dependences for the 10 MHz emission. 
It is obvious that the flattening is still present and thus cannot simply be 
caused by projection, but is an intrinsic feature of the emission. The emission 
pattern broadens as a whole (even in the direction perpendicular to the shower 
axis) as can be seen when comparing the back-projected patterns for a 45° 
inclined air shower depicted in Fig. 13 with that of a vertical shower shown in 
Fig. 3. 

The overall broadening of the emission pattern is due to the fact that the air 
shower maximum for inclined showers is much further away from the ground 
than for vertical showers. This effect was already predicted from geometri- 
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cal/qualitative arguments by j9|. It makes inclined air showers an especially 
interesting target for observation with radio techniques. 

The slight deviation of the 15° zenith angle curve from the trend seen in Fig. 
8 is explained by the shower's very small angle of only 5° to the 70° inclined 
geomagnetic field. The weakness of this deviation alone demonstrates that the 
dependence of the emission on the strength and orientation of the geomagnetic 
field is very slight — except regarding the polarization effects analyzed in Sec. 
4.2. Consequently, the same diagram for air showers coming from the north 
(not shown here) looks very similar. 

Fig. 10 shows the zenith angle dependence for the 55 MHz frequency compo- 
nent. The overall trend is the same as in the 10 MHz case, but the coherence 
losses cut off the emission pattern at a zenith angle dependent distance of a 
few hundred meters. At zenith angles ^30°, however, the coherence begins to 
hold up to large distances. This is confirmed by the spectra of a 45° inclined 
10^^ eV air shower shown in Fig. 11. The spectra are much flatter up to large 
distances when compared with the vertical case in Fig. 1. 

Inclined air showers thus not only offer significantly broader emission regions 
on the ground, but provide the advantage that the larger "footprint" even 
extends to significantly higher frequencies. 




1 00 200 300 400 500 600 700 800 900 1 000 
distance from shower centre [m] 

Fig. 8. Dependence of the 10 MHz east-west electric field component emit- 
ted by a 10^^ eV air shower coming from the south for different shower 
zenith angles as a function of distance to the north. Red/sohd: vertical shower, 
green/dashed: 15°, blue/dotted: 30°, violet/short dotted: 45°, turquois/dash-dotted: 
60°, black/double-dotted: 75° zenith angle. 
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Fig. 9. Same as Fig. 8 back-projected to the shower-based coordinate system (see 
text) . 
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Fig. 10. Same as Fig. 8 for the 55 MHz frequency component. 
4-2 Shower azimuth angle and polarization 



A very important trait of the radio emission is its predicted polarization, 
which is directly related to the shower azimuth angle. Knowledge of this de- 
pendence is imperative for the planning and interpretation of experimental 
measurements. 

Fig. 12 shows a comparison of the emission at 10 MHz from 10^^ eV air showers 
with 45° zenith angle as a function of azimuth angle. The total field strength 
pattern is elongated due to the projection effects arising at high shower zenith 
angles (cf. section 4.1). Taking out the pure projection effects leads to the 
patterns depicted in Fig. 13. The patterns are much more circular, but retain 
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Fig. 11. Spectra of the emission from a 10 eV air shower with 45° zenith angle 
at various distances to the north. From top to bottom: 20 m, 140 m, 260 m, 380 m 
and 500 m. 

a significant intrinsic ellipticity and asymmetry. 

The total field strength pattern of the emission (left column) simply rotates 
as a function of azimuth angle. (Deviations from a pure rotation are caused 
by the symmetry-breaking due to the magnetic field and shower axes — the 
emission pattern is no longer supposed to be truly symmetric.) In other words, 
no significant information associated with the geomagnetic field direction is 
present in the signal. As a direct consequence, it is not possible to verify 
the geomagnetic origin of the emission with an experiment measuring only 
the total field strength (or only one circular polarization component) of the 
emission. Furthermore, because the air showers arrive isotropically from all 
azimuthal directions, there will not be any azimuthal dependence of the mea- 
sured pulse amplitudes in statistical samples of measured total field strength 
pulses. It will therefore not be possible to confirm air showers as the source of 
measured radio pulses from statistics of total field strength data alone. In this 
case, independent information about the simultaneous arrival of cosmic rays, 
e.g. from particle detectors, would be necessary. 

The situation is different when one looks at the individual linear polarization 
components. (As demonstrated in section 3.4, the emission is polarized lin- 
early to a very high degree.) It is visible from the right three columns of Fig. 
12 that the field strength in the different polarization directions has a direct 
dependence on the geomagnetic field direction: the signal is linearly polarized 
mainly in the direction perpendicular to the air shower and magnetic field 
axes, at least in the central regions of high emission. The non-zero contribu- 
tions in the vertical polarization direction arise from the 70° inclination of 
the geomagnetic field. Fig. 14 illustrates the polarization characteristics of the 
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central emission region in a more intuitive way through indicators denoting 
the ratio of north-south to east-west polarization component over-plotted over 
the total field strength contours. 



Due to these polarization characteristics, experiments that measure the linear 
polarization characteristics of the emission can therefore directly verify the 
geomagnetic origin of the radio emission from cosmic ray air showers. 
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Fig. 12. Contour plots of the 10 MHz emission from a 10^^ eV air shower with 
45° zenith angle as a function of shower azimuth. Columns from left to right: 
total field strength, north-south, cast-west and vertical polarization component. 
Lines from top to bottom: 0°, 30°, 60° and 90° azimuth angle. Contour levels are 
0.25 fxV m-i MHz-i apart. 



4-3 Magnetic field 



As discussed in section 4.2, the magnetic field has an important influence on 
the polarization characteristics of the radio emission. The influence on the 
total electric field strength, however, is very weak. 

In Fig. 15 we compare the 10 MHz total field strength and polarization char- 
acteristics of vertical 10^^ eV air showers in four different magnetic field con- 
figurations: fields of 0.3 Gauss and 0.5 Gauss strength with horizontal and 
70° inclined geometry. A 0.3 Gauss horizontal magnetic field is present in the 
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Fig. 14. Contour plot of the total 10 MHz electric field component emitted by a 
10^^ eV 45° inclined air shower with over-plotted indicators denoting the ratio of 
east- west to north-south polarization. Left: azimuth of 0°, right: azimuth of 90°. 
Contour levels are 0.25 fiV m^^ MHz^^ apart. 

equatorial region, whereas a 0.5 Gauss ~70° inclined magnetic field is present 
in central Europe. 

The change from a horizontal magnetic field to a 70° inclined magnetic field 
introduces a number of effects. First, a small north-south asymmetry arises. 
Second, the overall emission level drops only very slightly — although the 
projected magnetic field that the vertical air shower sees drops by a factor of 
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cos~^(70°) ~ 3. This demonstrates the very weak effect of the magnetic field 
on the total emission field strength. The most prominent change is visible in 
the polarization characteristics along the east-west direction from the shower 
center. 

Increasing the field strength from 0.3 Gauss to 0.5 Gauss mainly boosts the 
fiux in the east-west direction from the shower center in case of a horizon- 
tal magnetic field. In case of a 70° inclined magnetic field, the changes are 
minimal. 
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Fig. 15. Contour plots of the total 10 MHz electric field emitted by a 10 eV 
vertical air shower with over-plotted indicators denoting the ratio of east-west to 
north-south polarization. Top: 0.3 Gauss magnetic field, bottom: 0.5 Gauss magnetic 
field, left: magnetic field horizontal, right: magnetic field 70° inclined. Contour levels 
are 0.25 jj^V m~^ MHz~^ apart. 

4-4 Primary particle energy 

Another important characteristic of the radio emission is its dependence on the 
primary particle energy. In Fig. 16 we present the dependence of the 10 MHz 
frequency component at various distances from the center of a vertical air 
shower as a function of primary particle energy. Although the depth of the 
shower maximum obviously depends on the primary particle energy, it is kept 
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constant in these calculations to assess the influence of the primary particle 
energy alone. (For a combined dependence see section 5.4). 



The scaling of the field strength with primary particle energy is approximately 
linear, following a power-law oc E^^^. For 55 MHz, the diagram is very similar 
until the curves again cut off at distances of a few hundred meters due to 
coherence losses (not shown here). The spectra do not change significantly in 
comparison with the 10^^ eV case except for an overall change in amplitude 
as demonstrated in Fig. 17 in comparison with Fig. 1. 

An approximately linear scaling of the emission with primary particle energy 
is expected for coherent emission. In the coherent regime, the field strength 
directly scales with the number of emitting particles. (The emitted power 
consequently scales as the number of particles squared.) Since the number 
of particles grows approximately linearly with primary particle energy in the 
parametrizations at the basis of our simulations, the linear scaling directly 
translates to a linear scaling of the field strength with primary particle energy. 
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Fig. 16. Scaling of the 10 MHz east-west electric field component emitted by a 
vertical air shower as a function of primary particle energy E^. From top to bottom: 
20 m, 100 m, 180 m, 300 m and 500 m from the shower center. The data follow a 
power-law oc E^'^^ . 



4-5 Depth of shower maximum 



The depth of the air shower maximum is directly related to the nature and 
energy of the primary particle [cf., e.g., l3|. Additionally, it is one of the pa- 
rameters that undergo strong fluctuations between individual showers with 
otherwise identical parameters. It is therefore interesting to evaluate the de- 
pendence of the emission on this parameter. 
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distances to the north. From top to bottom: 20 m, 140 m, 260 m, 380 m and 500 m. 

As can be seen in Fig. 18, there is a significant dependence of tlie emission on 
tlie deptli of the air shower maximum. The deeper penetrating the air shower, 
the steeper becomes the radial emission pattern. This is especially important 
for extremely high-energy air showers > 10^° eV, where the shower maximum 
can develop close to sea-level. The effect is the same as that visible in the 
zenith angle dependence (cf. section 4.1), where it is much more pronounced 
because the (spatial) distance of the shower maximum from the ground grows 
very rapidly with increasing zenith angle for a given value of X^s,x- 

The effect is very similar at 55 MHz, except for the cutoffs due to the loss of 
coherence at distances above a few hundred meters. 



5 Parametrizations for vertical showers 

Having analyzed the qualitative dependences of the radio emission on various 
air shower and observer parameters, constructing a parametrization of these 
dependences would be very useful. In a first step, we therefore quantify the 
dependences in a simple manner for vertical geometry. Afterwards, we gener- 
alize these dependences to an arbitrary geometry and piece together an overall 
parametrization taking into account the parameters simultaneously. 

As there is no direct error estimate for the underlying Monte Carlo results, 
we neither specify values nor make any error estimates for the derived 
fit parameters in the following sections. We specify our fit parameters with 
a large number of significant digits, knowing that the parameters are not 
determined with such high precision. Nevertheless, this allows an overall better 
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Fig. 18. Radial dependence of the 10 MHz component from a vertical 10^'^ eV 
air shower with various depths of the shower maximum Xmax- Red/solid: 
-'^max = 560 g cm~^, green/dashed: Xmax = 595 g cm~^, blue/dotted: 
-'^max = 631 g cm~^, violet/short dotted: Xmax = 665 g cm~^, turquois/dash-dotted: 
-^max = 700 g cm~^, black/double-dotted: Xmax = 735 g cm~^. 

representation of the Monte Carlo results with the fit functions. To verify the 
quality and estimate the deviation of our parametrization from the Monte 
Carlo data, we then make a direct comparison of our overall parametrization 
and the corresponding Monte Carlo results for a sample of test parameter sets 
in section 6.5. 



5.1 Radial dependence 



The radial dependence of the emission on distance r from the shower center 
can be fit with an exponential decay, 



E{r, 2Tciy) = Eq exp 



r 



This does not attempt to fit the behavior in the incoherent regime. Fig. 19 
shows the simulated 10 MHz and 55 MHz total field strength components as 
a function of distance to the north from the shower center with the associated 
exponential fits. For the 10 MHz component, we fit one exponential in the 
central 500 m and a second in the outer 500 m region. This increases the 
quality of both fits very significantly. For the 55 MHz component, we use 
only the values up to 380 m, as the emission becomes incoherent at larger 
distances. We do not take any asymmetry of the emission pattern (cf. Fig. 3) 
into account in this parametrization. 
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Table 1 

Parameters for the radial fits according to eq. (1) depicted in Fig. 19. 
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Fig. 19. Exponential radial dependence of the 10 MHz emission in the central 500 m 
(solid), 10 MHz emission in the outer 500 m (dotted) and 55 MHz emission (dashed) 
from a 10^^ eV vertical air shower with corresponding Monte Carlo simulated data. 

5.2 Spectral dependence 



The spectral dependence in the coherent regime can also be parametrized well 
with an exponential decay. The dependence in the incoherent regime at high 
frequencies is not well determined by the Monte Carlo simulations performed 
so far. It is, however, probably flatter than an exponential decay as demon- 
strated from the analytical calculations, where the functional form converges 
towards a power-law. Outside the valid frequency regime, the parametrization 
therefore is bound to underestimate the real flux. We flt the function 



E{r, 2t:u) = Eq exp 



[u - 10 MHz) 



(2) 



to spectra at various distances r from the shower center to the north. The 
parameter Eq in this case directly represents the total fleld strength at 10 MHz. 
Only data in the coherent regime is used for the fltting procedure. The data 
range used is indicated in Fig. 20 together with the resulting flt functions. The 
associated flt parameters are listed in Table 2. 
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0.7942 
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Table 2 

Parameters for the spectral fits according to eq. (2) depicted in Fig. 20. 




10 100 

[MHz] 

Fig. 20. Fit of exponential functions to spectra from a 10^^ eV vertical air shower 
at various distances to the north. From top to bottom: 20 m, 140 m, 260 m, 380 m 
and 500 m. Only the data in the coherent regime is used for the fit, as indicated by 
the colors. 

5.3 Polarization characteristics 



In the center region, the emission is almost purely linearly polarized in the 
direction perpendicular to the magnetic field and air shower axes (cf. Fig. 14), 
in good agreement with the analytic calculations. In particular, the electric 
field vector points in the direction 



sin 6 sin-i^B siny? 
cos 6 cos — cos (fi sin 6 sin -i^B 
cos '(9b sin^ sirup 



[cos 9 cos 't?B — cos (/p sin^ sin -(^b)^ + sin^ 6* sin^ ip 



(3) 
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where 9 denotes the shower zenith angle, is the shower azimuth angle with 
respect to the magnetic north and d-Q specifies the inclination angle (i.e., 
complement of the zenith angle) of the magnetic field. 

Multiplication of an E{r, 27nj, Ep) -value with the unit polarization vector 

E{9,ip,{}b) then directly yields the estimated north-south, east-west and ver- 
tical linear polarization components (in this order). 

The complex dependences at larger distances from the shower center cannot 
be easily parametrized at this stage. 



5.4 Combined Ep and X^ax dependence 



We have discussed the radio emission's dependence on the primary particle 
energy and the depth of the shower maximum separately in earlier sections. 
Here, we parametrize the combined dependence on primary particle energy 
and appropriately adjusted depth of shower maximum refiecting the deeper 
atmospheric penetration of higher energy air showers. We set X^ax to 500, 
560, 631, 700 and 770 g cm'^ for Ep values of 10^^ 10^^ 10^^ 10^^ and 
10^^ eV, respectively, see 12; 1^. The steepening of the radial dependence for 
increasing X^ax discussed in section 4.5 in this case leads to a radius-dependent 
steepening or fiattening of the energy dependence (originally oc E^'^^ as shown 
in section 4.4) in the central and outer regions, respectively. However, the 
combined dependence is still well- described by a power-law of the type 
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Fig. 21. Scaling of the 10 MHz electric field emitted by a vertical air shower as a 
function of primary particle energy Ep with appropriately changing depth of shower 
maximum Xmax- From top to bottom: 20 m, 100 m, 180 m, 300 m and 500 m to 
the north from the shower center. 
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%,27r^,Ep) =E, (j^y^^\ (4) 



,1017 eV 

The associated fit parameters are listed in Table 3 



r [m] 


£^0 [/^V m-i MHz-i] 


K 


valid Ep [eV] 


20 


10.18 


1.057 


10i5_ioi9 


100 


5.188 


1.004 


10l5_10l9 


180 


2.995 


0.965 


10l5_iol9 


300 


1.315 


0.907 


10l5_iol9 


500 


0.265 


0.808 


10l5_10l9 



Table 3 

Parameters for the combined primary particle energy and depth of shower maximum 
dependence according to eq. (4) depicted in Fig. 21. 



6 Parametrizations for arbitrary geometry 



We now generalize our parametrizations of the radio emission as a function of 
air shower and observer parameters to an arbitrary shower geometry. As our 
final result, we piece the individual parametrizations together to an overall 
parametrization incorporating all major parameters. 



6.1 Radial dependence 



The emission pattern becomes increasingly asymmetric with increasing zenith 
angle (cf. Fig. 12). As discussed in section 4.1, most of this asymmetry is 
caused by projection effects that can be taken into account by changing from 
a ground-based coordinate system (distance r from the shower center) to a 
shower-based coordinate system (perpendicular distance / from the shower 
axis), cf. Fig. 13. The remaining intrinsic asymmetries in the emission pattern 
we do not take into account in our parametrization. 

The back-projection from the ground-based to the shower-based coordinate 
system is given by 



l{r) = r^l- cos^ {ifio - (fi) sin^ {9), (5) 



when ip and 6 specify the shower azimuth and zenith angle and r and (po 
denote the observer distance from the shower center and the observer azimuth 



22 



angle (azimuth angles being measured with respect to the north). The 10 MHz 
back-projected radial dependence can then be well fit as 



E{1) = Ee exp 



(6) 



The fitting is performed as described in section 5.1 for each zenith angle 9 
individually. We restrain ourselves to the inner 500 m (back-projected) radius 
to increase the accuracy of the parametrization. We also base the fits on air 
showers coming from the north rather than the south to exclude the deviation 
for the 15° zenith angle case arising from the only 5° angle to the geomagnetic 
field in case of an air shower coming from the south. The resulting fits are 
depicted in Fig. 22, and the corresponding parameters are listed in table 4. 



^ 1 

0.1 

100 200 300 400 500 

distance from shower axis [m] 

Fig. 22. Radial dependence to the north for the 10 MHz emission from a 
10^^ eV air shower with corresponding exponential fits. Red/solid: vertical shower, 
green/dashed: 15°, blue/dotted: 30°, violet/short dotted: 45°, turquois/dash-dotted: 
60° zenith angle. 
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Table 4 

Parameters for the radial fits according to eq. (6) depicted in Fig. 22. 
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6.2 Spectral dependence 

Similarly as for the radial dependence, we now generalize our parametrization 
for the spectral dependence of the reference shower to an arbitrary shower 
geometry using the same exponential fits as adopted in section 5.2. Again, for 
each shower zenith angle individually, we fit a number of spectra at different 
distances / from the shower axis using the function 



We do not show the individual fits explicitly here. The resulting parameters 
Eg{l) and z/e(/) are tabulated in table 5. The dependence of z/g on the distance 
to the shower axis I can in turn be parametrized with an exponential function 

Ml) = e^'^'r (8) 



In fact, the parameter ag can be fixed to the same value for all cases of 6 
analyzed here at only minor loss of precision. The fits with ae fixed to a value 
of 47.96 MHz are shown in Fig. 23. The associated parameters for bg are listed 
in table 6. 

100 



^ 10 



1 

100 200 300 400 500 600 
distance from shower axis [m] 

Fig. 23. Parametrization of as a function of distance to the shower axis / ac- 
cording to eq. (8). Red/solid: vertical shower, green/dashed: 15°, blue/dotted: 30°, 
violet/short dotted: 45°, turquois/dash-dotted: 60° zenith angle. 
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Ee{l) exp 



[u - 10 MHz) 



(7) 
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Table 5 

Parameters for the spectral fits of 10^^ eV air showers with arbitrary geometry 
according to eq. (7). (The arbitrary-seeming values for I arc due to calculation of 
the air shower in the ground-based coordinate system with subsequent conversion 
to the shower-based coordinate system.) 
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Table 6 

Parameters for the parametrization of U0{1) according to eq. (8). 
6.3 Dependence of radial scale factor on Xmax 



To factor the influence of the (vertical equivalent) depth of shower maximum 
into the parametrization, we have to parametrize the flattening of the emis- 
sion's radial dependence with increasing Xmax- To achieve this, we fit the 
radial dependences calculated for vertical 10^^ eV showers with various values 
of ^max (shown in Fig. 18) with exponential functions in the central 500 m 
(not shown here). As can be seen in Fig. 18, the curves overlap at a distance 
of / = r ~ 200 m. Taking this point as a reference, the effect of changing Xmax 
can be reduced to a pure change of the slope, i.e., the scale factor of the 
exponential. We can then quantify the change of this scale factor by the ratio 

Oi{X^^) = —— -— (9) 

Zfl 631 g cm ^) 



of the scale factor for a given Xmax and the scale factor of our reference shower. 
Figure 24 shows o; as a function of Xmax and a fit of this dependence using a 
power-law 

/ Y \ -1-50519 

a (^max) = 1.00636 — '"^ J . (10) 

^ ' \631 g cm-2 J ^ ' 



6.4 Overall parametrization 



We now piece the individual parametrizations together to an overall formula. 
This imphes that the different effects are independent of each other and there- 
fore can be separated in an easy way, which need not be true in all cases. 
Nevertheless, such an overall parametrization can be a useful basis for com- 
parisons of experimental data with theoretical predictions as long as one keeps 
the limitations of the parametrization in mind. We therefore provide a formula 
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Fig. 24. Ratio a of the scale factor for a given X^ax to the scale factor of the 
reference shower as a function of X^ax and corresponding power-law fit. 

similar to the parametrization first given by jl^] and further enhanced by j^]- 
The associated sets of parameters are given for each zenith angle individually. 



At the heart of the parametrization is the radial dependence for arbitrary 
geometry as described in section 6.1. This we combine with the spectral de- 
pendence derived in section 6.2 and get 



— * 

E{r, v9o, z^) = f Ee exp 







exp 


[ le \ 



u/MRz - 10 



47.96 exp[-l{r,ip,)/be 



with the parameters as listed in table 7. The "fudge factor" / is introduced to 
allow deviations that might arise in the future when a more sophisticated 
distribution of particle momenta (i.e., pitch-angles) is taken into account. 
Throughout this work, we set / = 1. Additionally, we can factor in the primary 



e 


Eg [/iV MHz-i] 


Ig [ni] 


bg [m] 


0° 


12.33 


135.30 


219.41 


15° 


11.04 


152.80 


219.16 


30° 


8.33 


202.09 


254.23 


45° 


4.98 


339.71 


305.17 


60° 


2.53 


873.54 


590.03 



Table 7 

Parameters for the overall parametrization of the air shower emission according to 
eq. (12). 

particle energy dependence (cf. section 4.4) and the dependence of the radial 
scale factor on the (vertical shower equivalent) depth of the shower maximum 
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derived in section 6.3. For the latter, we have to take into account the change 
of the field strength in the shower center associated with the steepening of the 
scale factor Ig calculated from the reference point at / = 200 m. The resulting 
overall parametrization is then given by 



E(r, ipo, 1^, Ep, Xmax) Ee 

X exp 

X exp 



E^ 



0.96 



1017 eV. 
200m(Q;(X^a^)-l) + /(r, ipo) 

u/MRz - 10 
'47.96 exp [-Z(r, ^o)/fe0] 



;i2) 



with a(Xinax) defined in eq. (10), l{r,ipo) given by eq. (5), the values for Eg, 
le and bg taken from table 7, and / = 1. To calculate the field strength for the 
individual linear polarization components, one can then multiply the result of 
eq. (12) with the unit polarization vector given in eq. (3). 



6. 5 Quality and validity of the overall parametrization 



To verify the quality of our parametrization, we take a sample of observer 
parameters testing the different regimes of the parametrization and compare 
the result of eq. (12) with the result from the Monte Carlo simulation. As can 
be seen in table 8, the deviations are acceptably small: most of the time the 
error is below 10%, and only in cases where the parametrization is expected to 
degrade (e.g. in the east-west direction of heavily inclined air showers, where 
the intrinsic asymmetries in the emission pattern become relevant) it grows 
beyond 20%. 

One must of course be careful not to leave the parameter regimes for which the 
parametrization was created. Specifically, the back-projected radial distance 
/ was limited to 500 m in the underlying radial fits. The frequency limits 
for the spectral fits can be estimated from table 5 as a function of I. As 
explained earlier, the parametrization is bound to severely underestimate the 
fiux at higher frequencies. The polarization characteristics given by eq. (3) 
are only valid for the central region as illustrated by Fig. 14. At significant 
zenith angles, the intrinsic asymmetries of the emission pattern which were 
not taken into account in the parametrization lead to a growing deviation 
from the Monte Carlo results. Finally, special caution should be used when 
changing the depth of shower maximum X^ax for significantly inclined air 
showers, as the underlying parametrization of Q;(Xmax) was derived for vertical 
air showers only and the projection effects associated with inclined showers 
greatly enhance the depth of shower-maximum effects. 
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Table 8 



Quality check of the overall parametrization given by eq. (12). 
6.6 Comparison with AUan-parametrization 



The authors of [lj| provided a parametrization of their experimental data 
which has in turn been generalized by to the form 



= 13 /iV m-i MHz-i 
X exp 



—r 



1017 eV 

-1 



50 MHz 



sma cos I 



sin 45° cos 30"^ 



(13) 



where a denotes the angle between shower axis and magnetic field, 9 is the 
shower zenith angle and tq is a scale factor of about 110 m. To compare this 
with our results, we have to convert the experimentally motivated values to 
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our theoretically derived E{uj) values. For the conversion we use the relation 
E{R,uj) ^ 6.4 E{R,uj) (14) 



'128 



TT 



as derived in |5|]. The conversion results in 



E^ 



sin a cos 6* 



10^7 eVj V sin 45° cos 30° 



X exp 



— r 



50 MHz 



(15) 



Our parametrization yields a value of 3.6 fiV MHz~^ for the 50 MHz 
emission in the center of a 30° zenith angle 10^^ eV air shower, which is not 
far off the Allan- value of 2 /iV MHz~^. (For comparison, the estimated 
theoretical limit for a 3a detection with 1/10/100 LOPES antenna(s) corre- 
sponds to ~ 0.4/0.15/0.05 /xV m~^ MHz~^, see jsf.) The radial dependence 
in both cases is given by an exponential decay, and the resulting scale factor 
in our parametrization indeed corresponds to ~110 m for the aforementioned 
set of parameters. The linear scaling with primary particle energy is identical. 
Apart from these similarities, there are of course some differences: the fre- 
quency dependence in the Allan-formula is specified as u'^. This is obviously 
very different from the exponential decay in our parametrization. However, 
the extrapolated frequency dependence for the Allan-formula rests on rather 
sparse and uncertain data, some of which lie in the incoherent regime not 
included in our parametrization. The dependence on shower zenith angle is 
much more complex than a simple cos(^) trend in our parametrization and 
is therefore difficult to compare with the Allan-formula. Furthermore, we do 
not predict any significant dependence of the total field strength on the angle 
between shower axis and magnetic field (denoted a in the Allan-formula). The 
emission in an individual linear polarization component, on the other hand, 
varies with the angle between shower axis and magnetic field axis (cf. eq. (3)). 

Overall, our parametrization shows many similarities to the historic Allan- 
parametrization. The discrepancies, e.g. regarding the frequency dependence, 
are significant, but considering the sparse experimental data on which the 
Allan-formula is founded, these discrepancies should not be over-interpreted. 
In this context one should also remember that later experiments measured 
significantly lower values for e^, a discrepancy that is yet unsolved, but most 
probably due to calibration issues. 

On one hand, we therefore urgently need new, reliable, well-calibrated exper- 
imental data. On the other hand, the modeling efforts have to continue. In 
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particular, we will improve our simulation by basing it on a more realistic 
air shower model as given by the CORSIKA code. This will automatically 
resolve the two major shortcomings of our current model: an unrealistic parti- 
cle pitch-angle distribution and the assumption of a totally homogeneous air 
shower development. Taking into account these effects most hkely will redis- 
tribute flux in such a way that the emission levels in the center region become 
smaller, making them more consistent with the historical data. 



7 Discussion 

Our analysis of the Monte Carlo simulation results with regard to the under- 
lying air shower parameters for the first time establishes a number of experi- 
mentally relevant features of the radio emission. 

One major result is the predicted polarization characteristics of the emis- 
sion generated in the geomagnetic emission scenario. With this knowledge, 
polarization-sensitive experiments should be able to directly verify that a ma- 
jor part of the radio emission indeed stems from the geomagnetic mechanism. 
The overall weak intrinsic asymmetries in the emission pattern (except for 
those associated with projection effects) make experimental setups measuring 
only the total field strength or one circular polarization component seem less 
desirable. 

Another important insight is provided by the effects arising in air showers with 
large zenith angles. The intrinsic broadening of the emission pattern associated 
with the increasing distance of the air shower maximum in combination with 
projection effects and a flattening of the spectral dependence makes highly 
inclined air showers an especially interesting target for detection with radio 
techniques. 

Other useful results of our simulations are the predicted wavefront curvature, 
the expected quasi-linear scaling of the field strengths in the coherent regime, 
the very weak dependence of the total field strength on the specific geomag- 
netic field geometry and strength, and the changes to the radial dependence 
as a function of changing depth of shower maximum. 

The successful incorporation of a significant number of air shower and observer 
parameters into a single parametrization as achieved in this work demonstrates 
that the emission is overall "well-behaved" and yields a well-interpretable 
signal. As a solid estimate, the parametrization can be a useful tool for the 
interpretation of experimental data and the planning of experimental setups. 
Additionally, there are remarkable similarities between the parametrization of 
our Monte Carlo results and the historic Allan-formula, although significant 
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uncertainty regarding the absolute calibration of the historical data remains. 

Our current results represent the most sophisticated simulation of radio emis- 
sion from cosmic ray air showers carried out to date. In the future, we will 
switch from analytic parametrizations of air shower characteristics as the basis 
of our simulations to a full-fledged air shower model based on the CORSIKA 
simulation code, further improving the modeling accuracy. The two main draw- 
backs of the current model, an over-simplified particle pitch-angle distribution 
and the disregard of inhomogeneities in the air shower development, will then 
automatically be resolved. 



8 Conclusions 



This work presents the cumulative result of our effort at a realistic modeling 
of the radio emission from cosmic ray air showers in the scheme of coherent 
geosynchrotron radiation. For the first time, we now have a solid understand- 
ing and a quantitative description of the important emission characteristics 
and their dependences on important air shower and observer parameters. In 
particular, the emission pattern, spectral dependence, polarization character- 
istics, primary particle energy dependence, magnetic field dependence and 
the dependence on air shower geometry and depth of shower maximum for 
geosynchrotron emission are now theoretically determined. This information 
is imperative for the interpretation and planning of concrete experiments. In 
the near future, on the other hand, experiments such as LOPES will provide 
well-calibrated, reliable data that will allow a direct comparison with our the- 
oretical predictions for the first time. A direct verification of the geomagnetic 
emission mechanism will then be possible. 

Having reached this milestone in the modeling of radio emission from cosmic 
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ray air showers, the next step will be to interface our code with CORSIKA 
delivering a full-fledged realistic Monte Carlo simulation of radio emission from 
cosmic ray air showers. Additionally, we wish to incorporate further possible 
emission mechanisms such as Askaryan-type Cerenkov radiation in our model. 
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